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ABSTRACT 

The X-ray variations of hard state black hole X-ray binaries above 2 keV show 'hard 
lags', in that the variations at harder energies follow variations at softer energies, with 
a time-lag r depending on frequency i/ approximately as r oc u~^-'^ . Several models 
have so far been proposed to explain this time delay, including fluctuations propagating 
through an accretion flow, spectral variations during coronal flares, Comptonisation 
in the extended hot corona or a jet, or time-delays due to large-scale reflection from 
the accretion disc. In principle these models can be used to predict the shape of the 
energy spectrum as well as the frequency-dependence of the time-lags, through the 
construction of energy-dependent response functions which map the emission as a 
function of time-delay in the system. Here we use this approach to test a simple re- 
flection model for the frequency-dependent lags seen in the hard state of GX 339-4, 
by simultaneously fitting the model to the frequency-dependent lags and energy spec- 
trum measured by XMM-Newton in 2004 and 2009. Our model cannot simultaneously 
fit both the lag and spectral data, since the relatively large lags require an extremely 
fiarcd disc which subtends a large solid angle to the continuum at large radii, in dis- 
agreement with the observed Fe Ka emission. Therefore, we consider it more likely 
that the lags > 2 keV are caused by propagation effects in the accretion fiow, possibly 
related to the accretion disc fluctuations which have been observed previously. 

Key words: accretion, accretion discs — black hole physics — stars: individual (GX 
339-4) — X-rays: binaries 



1 INTRODUCTION 

Black hole X-ray binaries (BHXRBs) have extensively been 
studied both in terms of their spectra (where a soft, multi- 
colour disc black-body and a harder power-law component 
play the main role) and in terms of their variability, which 
can be straightforwardly quantified using the common ap- 
proach of computing the Power Spectral Density (PSD) of 
their signal. In combination with the time-averaged X-ray 
spectrum, this approach has enabled the variability to be 
studied as a function of spectral state, so that the hard 
power-law-dominated states can be associated with large 
variability amplitudes (tens of per cent fractional rms) and 
band-limited PSD shapes, while the soft disc-dominated 
states show much weaker variability and broadband power- 
law like PSD shapes. 

Although they provide information on the global con- 
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nection between the emitting components and the variability 
process, the PSD and mean energy spectra alone carry no in- 
formation about the complex pattern of intcrlinkagc between 
diverse regions in the accreting system that lead to transfer 
of variability from one physical component to another over a 
range of time-scales. In order to understand the causal con- 
nection between emission processes that show different rel- 
ative strengths in two difi'orcnt energy bands, the time-lags 
between these two bands can be extracted as a function of 
Fourier frequency. To date, frequency-dependent lags have 
been studied mostly in the hard state of BHXRBs. There 
the lags (r) are 'hard', in the sense that variations in harder 
bands lag behind variations in softer bands and depend on 
frequency as t oc v'"'^ , albeit with some sharper 'steps' in 
the lag-frequency relation ([Miyamoto et al. 19881 ICui et af] 
1997; Nowak et al. 1999). The general form of these lags has 
been explained by a variety of models, including Compton- 
isation in the extended corona (Kazanas et al. 1997) or a 
jet | [Reig et al.|2003(|Kylafis et al.|2008| ), spectral variability 
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during coronal flares (Poutancri & Fabianj[1999|, accretion 
fluctuations propagating through a power-law emitting re- 
gion (e.g. a corona) with an energy- dependent radial emis- 
sivity proflle ( |Kotov et al.||2001i [Arevalo fc Uttley[|2006|), 



or light-travel times to an extended reflecting region | |Kotov 
et al.|[200T Poutanen 2002). Models invoking Comptonisa- 



tion in the extended regions predict broader auto-correlation 
function for photons at higher energies, which suffer more 
scattering, which is opposite to what is observed ( Maccaronej 
|et al.|2000||Poutaneii|2001||. 

Recent measurements of the lag (Uttley et al. 2011) be- 
tween variations of the accretion disc blackbody emission 
and the power-law component strongly indicate that the 
variations are driven by fluctuations propagating through 
the disc (|Lyubarskii||1997| to the power-law emitting re- 
gion, however this model cannot simply explain the hard 
lags between bands where the power-law dominates. 

Models to generate frequency-dependent lags work by 
determining the response to an input signal of an emitting 
region which produces a hard spectral component, e.g. a 
compact coronal region sandwiching the disc, the upscatter- 
ing region in the jet or an extended reflector. The input sig- 
nal may be changes in the geometry, a fluctuation in accre- 
tion rate, seed photon illumination or illuminating primaxy 
continuum. The hard emitting region can only respond after 
a delay which, broadly speaking, is set by the time taken for 
the signal to propagate to and across the region. The delay 
time is determined by the signal speed (viscous time-scale 
or light-travel time) and the size scale of the hard emitting 
region. If the delay is large compared to the variability time- 
scale, the variations of the hard emitting region are smeared 
out and so the amplitude of variations of the lagging compo- 
nent is reduced. Thus the observed drop in lags with Fourier 
frequency can be reproduced. 

In principle, models for the emitting regions which can 
reproduce the lags should also be able to reproduce the en- 
ergy spectrum. Therefore, the combination of lag informa- 
tion with information on the X-ray spectral shape should 
provide rrmch greater constraints on models for the omitting 
regions than either the commonly used spectral-fltting meth- 
ods or the rarely-attempted fits to timing data. As a good 
proof-of-principle, the simplest models to attempt this joint 
lag and spectral-fitting are those where the lags are produced 
by reflection, i.e. so-called 'reverberation' lags. Evidence for 
small (few to tens of i?g light-crossing time) reverberation 
lags generated by reflection close to the black hole has been 
seen in Active Galactic Nuclei (Fabian & Ross 2010, Zoghbi 
et al. 2010: de Ma rco et al. 2011^ .Emmanoulo poulos et al. 
[2011) andBHXRBs (Uttl ey et al.|201ip , however the much 
larger lags we consider here require a larger scale reflec- 
tor, perhaps from a warped or flaring outer-disc, which was 
originally considered by Poutanen (2002). Reflection of hard 
photons from a central corona off an accretion disc should 
produce both a reflection signature in the spectrum and a 
frequency-dependent time-lag, the shapes of which should 
depend on the exact geometry of the reflector as well as 
on the location of the corona itself (for a central, point-like 
corona, this would be the height above the disc). It has to 
be stressed that only non-flat disc geometries are expected 



XMM-Newton provides a time resolution of 5.965 ms (frame 
readout time) in Timing mode with a 99.7% livetime and an 
energy resolution ~ 125 eV at 6 keV that combined, offer 
great potential for developing models that take into account 
energy and variability information together down to mil- 
lisecond time-scales. Our aim in the present paper is to fit 
the spectra and time lags from an EPIC-pn Timing mode 
observation of a BHXRB simultaneously to discern whether 
reflection is the main driver of the observed hard-to-medium 
lags. In order to achieve this, we have developed reflags, 
a reflection model that assumes a flared accretion disc act- 
ing as reflector (assuming the constant density ionised disc 
reflection spectrum of jBal lantyne et al. 2001), and is able to 
output for a given geometry either the resulting spectrum, 
or the expected lags as a function of frequency. We use this 
model to fit simultaneously the spectra and time-lags of the 
low-mass BHXRB GX 339-4 in the XMM-Newton obser- 
vations of the hard state obtained in 2004 and 2009, and 
so determine whether the refiection model can explain the 
hard lags observed in these observations while also remain- 
ing consistent with the X-ray spectrum. In principle, this 
combined approach could yield much greater constraints on 
the outer disc geometry than can be obtained with spectral- 
fitting alone, since low velocities of the outer regions of the 
disc cannot be resolved with CCD detectors, although these 
regions should contribute significantly to the lags if they 
subtend a large solid angle as seen from the source. 

We describe our model in Section |2] and the data reduc- 
tion and extraction of spectra and time-lags in Section 3 In 
Section l4lwe discuss the results from fitting the energy spec- 
tra and frequency-dependent time-lags of GX 339 4 together 
with our model, and also compare the expected optical/UV 
reprocessing signature of the inferred geometries with data 
from the Swift satellite. We discuss our results and the wider 
implications of our combined spectral-timing model-fitting 
approach in Section |5j 



2 A FLARED ACCRETION DISC MODEL 
2.1 Model parameters 

We consider a simple refiection model where the hard lags 

are produced by the light travel times from a variable power- 
law continuum source to the surface of a fiared accretion 
disc which absorbs and reprocesses the incident radiation or 



scatters it producing a hard refiection spectrum | George & 
|Fabianl|1991l). Following the geometry described by |Pouta- 
nen (2002), we have written reflags, an xspec model that 
can describe both the lags as a function of frequency and 
the mean spectrum, allowing simultaneous fitting by tying 
together the parameters from the spectral and lag fits. We 
assume that the accretion disc is axially symmetric but has 
a height of the disc surface above the mid-plane which de- 
pends on radius r as a power-law 



z{r) = Hout 



(1) 



to contribute a significant lag at low frequencies I jPoutanen 



20021. 



The CCD technology of the EPIC-pn camera onboard 



where i^out is the height of the disc surface at the outermost 
radius Rout and 7 is the 'fiaring index'. 

We place a point-like source of Comptonised photons at 
a certain height i?src above the disc, located in the axis of 
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symmetry of the system, and assume that the disc is trun- 
cated at some inner radius -Rin. Because the lags are propor- 
tional to the distances, while the spectral distortions due to 
rapid rotations are a function of radii in units of gravita- 
tional radius Rg = GM/(? , we have to specify the black 
hole mass, which we assume to be Mbh = lOM©. 

The amount of reflection of Comptonised photons that 
is expected to come from each region of the disc strongly 
depends on the value 7, and the only case where the contri- 
bution of outer radii to spectra and lags can be significant 
is for concave (7 > 1) geometries (Poutanon 2002). 

The model also depends on the two parameters that 
describe the spectral shape: the incident power-law pho- 
ton index F and the ionisation parameter 5. Due to com- 
putational constraints in the model evaluation, the latter 
is approximated to be constant throughout the disc. For 
computational purposes, the disc is also divided into 100 
equally spaced azimuthal angles and 300 logarithmically 
spaced radii, to produce of a total of 30 000 cells. 

The total energy-dependent power-law plus reflection 
luminosity per unit solid angle at a time t that is emitted 
by a system whose geometry is described by the parameters 
above (represented by the set a) and seen at an inclination 
angle i with respect to the axis of symmetry, equals: 



(2) 



where the sum is performed over the index I that repre- 
sents a single disc cell, Lpi^{E,t) is the luminosity per unit 
solid angle produced by power-law emission and a{E) is the 
albedo function. 

The factor ita^i{E,i) contains the several projection 
terms and solid angle corrections required for the l-th cell 
for the geometry described by a (see [Poutanen 2002), as 
well as relativistic Doppler and gravitational redshift cor- 
rections (following the same simplified approach as the 
DISKLINE model, Fabian et al. 1989), and is normalised so 
that iia,i{E,i) = 1. For a given geometry a and incli- 
nation angle i, fiefi then equals the solid angle subtended 
by the disc, corrected for the source inclination as seen by 
the observer. Finally, n represents the time delay between 
observed direct power-law emission and refiected emission 
coming from the l-th cell. 

2.2 System response and time-lags 

For a direct power-law emission pulse (a delta-function) 

equation ( 2 ) can be rewritten as: 

Ltot{E,t\a,i) = Lpl{E) 

.a 



{S{t) + a{E)'-^ Yl ^o.AE,i)6{t - r,)) , 



(3) 

and can be factoriscd as 

Lt^t{E,t\a,i) = Lvi^{E)T{E,t\a,i), (4) 

where 

T{E,t\a,i)=5(t) + a{E)^^Y Ko.,i{E,i)5{t-Ti). (5) 



The system response function T[E,t\a,i) contains the 
energy-dependent time redistribution of an input power-law 
luminosity as a function of time in a geometry a, as observed 
at an inclination i. Provided that the energy spectrum is av- 
eraged over a time significantly longer than the largest r;, 
Tmax (which iu this case is of order 10 s), the spectrum is 
given simply given by: const, x /J^"""' T(E, t'|a, i)At' . By tak- 
ing the Fourier transform of equation Q: 



Uot{E,v,i) = T{Lv^{E)T{E,t\a,i)\ 
= Lvi.(E)f{E,v\a,i). 



(6) 



In general, the time-lags between two light curves s(t) 
and h(€) in the soft and hard bands can be calculated by 
computing their Fourier transforms S(y) — J^[s{t)] and 
H{i') = T[h{t)] and forming the complex-valued quantity 
C(i/) = S*{i')H{i'), called the cross-spectrum (where the 
asterisk denotes complex conjugation). Its argumen t is the 
phase lag or phase difference between S (u) and Hiv) I Nowak 
et al.|1999| |: 

<t>iu) = arg [Ciu)] = arg [(S* (i^)^(i^))] (7) 



and therefore 



(8) 



equals the frequency-dependent time-lag. 

The transfer function T{E, u\a, i) is equivalent to the 
response function T{E,t\a,i) in the Fourier-frequency do- 
main. It is then possible to take the expression above to com- 
pute the cross-spectrum that will give the time-lags caused 
by reflection between two broad energy bands s and h and 
obtain the lags as above (the notation has been simplified): 



C{u) = Lp-^^{E,)Lp^{E^)f'{E,,v)f{E^,v), 



r{y) = 



arg[T*(£;„z/)T(Sh,i/)] 



27rz/ 



(9) 



(10) 



3 DATA REDUCTION 

3.1 Extraction of spectra, event files and 
instrumental response files 

GX 339-4 was observed using the XMM-Newton X-ray 
sateUite on 2004 March 16 and 2009 March 26 (Obslds 
0204730201 and 0605610201, respectively). In the present 
work, we will concentrate on the data taken using the EPIC- 
pn camera in Timing Mode, which allows for a fast read- 
out speed, as opposed to imaging modes. This is done by 
collapsing all the positional information into one dimension 
and shifting the electrons towards the readout nodes, one 
macropixcl row at a time, thus allowing a fast readout with 
a frame time of 5.965 ms (Kustcr et al. 2002). 

The raw Observation Data Files (ODFs) were obtained 
from the XMM-Newton Science Archive (XSA) and reduced 
using the XMM-Newton Science Analysis System version 
10.0.0 tool epproc using the most recent Current Calibra- 
tion Files (CCFs). 

With the aid of evselect, the source event informa- 
tion was then filtered by column (RAWX in [31:45]), and 
pattern (PATTERN <= 4, corresponding to only single and 
double events), and only time intervals with a low, quiet 
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background (where PI in [10000 : 12000] ) were selected for 
the subsequent data analysis. The total exposure times are 
~ 164 ks and ~ 39 ks for the 2004 and 2009 observations, 
respectively. The background spectrum was obtained from 
Obsid 0085680501 between columns 10 - 18 when the source 
was fainter so that the background is not contaminated by 
the source (see Done & Diaz Trigo (2010) for further details), 
rmf gen and arf gen were used to obtain the instrumental re- 
sponse files. To account for the effects of systematic uncer- 
tainties in the instrumentation, an additional 1% was added 
to the error bars. 

In the present work, uncertainties in the estimation of 
the parameters are quoted at the 90% confidence level for 
one parameter of interest. 

3.2 Extraction of the time-lags 

When measuring the lags for real, noisy data one needs to 
first take the average of the cross-spectrum over many inde- 
pendent light curve segments (and also adjacent frequency 
bins). This is because observational noise adds a compo- 
nent to the cross-spectrum which has a phase randomly 
drawn from a uniform distribution (see e.g. Nowak et al. 
|l999). By averaging over many independent measures of the 
cross-spectrum, the contribution of these noise components 
can be largely cancelled out (the residual error provides the 
uncertainty in the lag). Therefore, the time-lags t{v) were 
extracted from the argument of the cross spectrum aver- 
aged over many segments of the light curve as the exposure 
time and the dropouts due to telemetry saturation permit. 
Throughout this work we will take as the soft band the in- 
terval 2.0-3.5 keV and as a hard band the interval 4.0-10.0 
keV. 

We choose to have 8192 bins/segment for a total dura- 
tion of ~ 24.43 seconds per segment, giving 6690 segments 
for 2004 and 1582 segments for 2009. This will also constrain 
the frequency ranges over which lags can be obtained. These 
frequencies have also been rebinned logarithmically with a 
step of A In = 0.14 in order to improve the signal-to-noise 
per frequency bin. 

The uncertainties on the lag measurements follow from 
the technique in Nowak et al. (1999) (see also Bendat & 
Piersol 2010). These scale as \/NM where A'^ is the number 
of segments used to average the cross-spectrum and M is 
the number of frequencies averaged per bin. The lags are 
plotted on Fig. 1, Given the fact that, within the errors, the 
time-lag dependence on Fourier frequency is consistent be- 
tween the two observations, we will sometimes assume that 
the time-lags for the two observations are equivalent and do 
not vary between the two observations as their shapes can- 
not be distinguished within the errors. This 'substitution' 
should in principle give us tighter constraints for our study 
in Section 4 



4 RESULTS 

We test the model described in Section 2 using the 2004 
and 2009 observations of GX 339-4, in order to understand 
the validity of the model and infer the geometrical param- 
eters of the reflecting disc that is required to explain both 
the spectrum and the lags. To do this, the same model is 




Frequency [Hz] 



Figure 1. Time-lags between 4.0-10.0 and 2.0—3.5 keV energy 
bands as a function of frequency for the 2004 (black) and 2009 
(red) observations. 




1.05 



0.95 



Energy [keV] 



Figure 2. Ratio to a powerlaw after fitting the 2.0—5.0 keV and 
7.0-10.0 keV continuum (black: 2004, red: 2009). The presence of 
the iron line is clear. 



concurrently fitted to the spectrum as well as the lag ver- 
sus frequency data. Using this approach, it is possible to 
discern the importance of reflection to explain the observed 
lags, and whether or not an extra variability component to 
produce the lags is required. 

In Fig. 2 we plot the ratio of the 2004 and 2009 spectra 
to a power-law fitted to the 2.0-10.0 keV interval excluding 
the region 5.0-7.0 keV. The iron line shape appears clearly 
different, with the 2004 dataset showing a more broadened 
profile that is skewed towards higher energies, while the 2009 
line appears narrower and peaked at about 6.4 keV, the value 
expected for neutral or weakly-ionised emission. 

We perform fits using the 2004 and 2009 spectra in the 
3.0-10.0 keV band as well as their respective lags (we name 
these model fits A and B, respectively). These energy in- 
tervals are chosen in order to avoid contamination from the 
band dominated by the disc blackbody emission. Due to the 
small uncertainties on the spectral data, the weak steepen- 
ing in spectral shape seen below 3 keV can skew the fit to 
the spectrum, hence we cut off the spectral fit at 3 keV. 
However, as shown by Uttley et al. (2011), the lags are not 
significantly affected by the disc at energies down to 2 keV, 
hence we include photons down to this energy in the lag 
determination, in order to increase signal-to-noise. 

The best-fitting parameters for fits A and B are shown 
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Parameter Obs. A Obs. B 



Rin/Rg 


< 130 


< 220 


Rout 1 Rg 


(16.8l^;^) X 10^ 


(20l^°) X 103 


J^out/ Ji-out 


^ u.oo 




Hgrc/ Rg 


> 380 


> 330 


7 


2.18l»i« 


> 1.69 


i 


39.7t°;t 


29.8I4I 


V 


-^-0.002 


1.4b8_(, 


log 5 


Z.UO_Q 02 


2 nq+o oa 
^■'^"^-0.01 


n<=«7(27r) 


1.13 


1.28 


(spectrum) 


1405 


1305 


(time-lags) 


152 


32 




1557 / 1417 


1337 / 1416 



Table 1. Spectral and lag fit parameters using the model REFLAGS 
in the band 3-10 kcV. Hereinafter, we quote the total as well as 
the contribution of the that corresponds to both the spectrum 
and lags, separately. 



in Table [Ij The corresponding model comparisons to the 
data, including the data-to-model ratios are plotted in Figs[3| 
and|4| These initial fits suggest that the model can represent 
the 2009 (B) data better. The upper limit on the disc inner 
radius is below 220 Kg in both cases. This value does not 
give us any information on whether the accretion disc in 
GX 339-4 is truncated or not (for discussions about disc 
truncation, see e.g. Tomsick ct al. 2009). 

As for the values of the outer radius, its value could 
reach up to ~ 8 x 10'* Rg in B, rrmch larger than what is 
found in A, indicating a larger contribution to the narrow 
iron line in this case. This is in agreement with the stronger 
core of the line in B as seen in the ratio plot (Fig. 2 ) as well as 
the difference between the solid angles. An extremely high 
value for the source height i?src and a high value for the 
tlout I Rout ratio suggest that the fits are being driven by the 
lags, whose amplitude is strongly dependent on the distances 
to the furthest regions of the reflector, even if their contribu- 
tion is small. The ionisation parameter remains consistent 
between the observations, although a visual inspection of the 
residuals and a clear difference in are indicative of a disc 
that can be described with a line of ionised iron, and whose 
contribution would likely come from the inner disc regions 
of the accretion disc where the ionisation could plausibly be 
larger. The blue wing of a relativistically broadened disc line 
could conceivably also contribute to the residuals which we 
fit with a narrow line, however since the inner disc radius is 
not strongly constrained by the lags, the model is already 
relatively free to fit this feature with relativistic emission, 
by being driven by the spectral data alone. The fact that 
it does not fit these residuals suggests that a more complex 
ionisation structure of the disc is likely (also see Wilkinson 

| 20TT| ». 

However, as seen in the right panel of Fig. the model 
is clearly under-predicting the lags that are observed by a 
factor of ~ 2, mainly at low frequencies that would corre- 
spond to a light-crossing time expected from distant refiec- 
tion. Given that our focus is on the lags and our aim is to 
understand whether they are compatible with being caused 
by reflection, we test whether the lags are being constrained 
by the line shape by adding an extra Gaussian component of 



Parameter Obs. A -|- extra line Obs. B' 



-Rin/-Rg 


< 90 


< 120 


■Rout jRg 


(8.815:6) X 10^ 


(ll.OlJ-6) X 103 


Hout 1 Rout 


0.84 ± 0.04 


U.90_(, ^4 


Hsvcl Rg 


4301^0 


> 380 


7 


2 29+0.21 
^■^^-0.33 


9 fi9+0-36 
2-b2_p 43 


i 


48 9+^^" 

-O.fi 


-^^•■'-0.8 


r 


1.479 ± 0.003 


1.46 ±0.05 


log? 


1 50+°l2 
1.0U_Q 24 


1 fiC+O.OS 
J--OO_0.17 


Ec [keV] 


6.92 ±0.02 




o-c [keV] 


Q 1 7-H0.03 
U-J-'_o.02 




n<=ff/(27r) 


0.78 


1.03 


X^ (spectrum) 


1132 


1335 


X^ (time-lags) 


107 


77 


xVdof 


1239 / 1414 


1412 / 1417 



Table 2. Spectral and lag fit parameters using the model REFLAGS 
plus an extra Gaussian fitted to 2004 spectrum and lags (obs. 
A), and mixing the 2009 observation spectrum with the 2004 
observation lags (obs. B'). 



variable width to the spectral model for the 2004 data only, 
that corresponds to the dataset with higher signal-to-noiso. 
In this way, we improve the fit 'artificially' so that the model 
can find the best-fitting contribution to the narrow 6.4 keV 
core which might originate at larger radii, and thus we allow 
the model more freedom to produce larger lags which better 
fit the data. 

The best-fitting parameters after this procedure can be 
found in the first column of Table [2j and an improvement in 
is clear. The centroid energy for the additional emission 
line does not correspond to any ionised iron fluorescence 
transition, hence its value may come from a weighted mean 
of lines with a complex ionisation. 

The decrease in the reflags solid angle is compatible 
with the addition of the extra line, i.e. the main model does 
not need to fit it anymore. Despite the attempt to artiflcially 
improve the spectral fit to allow the model to better fit the 
large scale reflection, the lags from the model are still too 
low (Fig.[5|. 

Given the fact that the lags in 2004 and 2009 are very 
similar despite a factor of 2 smaller error bars in 2004, we 
also try to flt the 2009 spectrum swapping 2009 and 2004 
lags (B', hereinafter). This worsens the fit with respect to 
fit B (Fig. 6 , second column in Table 2 ) , as the smaller error 
bars in the 2004 lags push the geometrical parameters to 
more extreme values, resulting in an increase of the resid- 
uals around the iron line and a slightly larger width of the 
line due to the increase of the flaring parameter accompa- 
nied by a decrease of the outer radius. These changes also 
increases the frequency at which the lag-frequency depen- 
dence becomes constant, which corresponds to the longest 
time-scale of variations in the response function). 

The disc parameters in all these fits lead to a geometry 
where the source height, disc flaring as well as the Hout I Rout 
ratio at the outer radius are extreme. This is due to the fact 
that the lags are driving the fits towards increased outer 
reflection to increase their amplitude, however the spectral 
model is not necessarily sensitive to small X-ray flux con- 
tributions from outer radii or the change in iron line shape 
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Energy [kcV] Frequency [Hz] 

Figure 3. Obs. A (2004) spectrum (left) and lags vs frequency (right). The thick solid line represents the model for the best-fitting 
parameters (see Table 1, left column). In the spectrum, the thin solid line represents the direct continuum power-law. The lower panels 
show the ratio data/model. 




Energy [keV] Frequency [Hz] 

Figure 4. Obs. B (2009) spectrum (left) and lags vs frequency (right). The thick solid line represent the models REFLAGS for the 
best-fitting parameters (see Table 1, right column). 
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Figure 5. Obs. A (2004) spectrum (left) and lags vs frequency (right). The thick solid lines represent the model REFLAGS (including a 
Gaussian component in the case of the spectrum) for the best-fitting parameters (see Table [2j left column). 
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Figure 6. Obs. B' (2009 spectrum, 2004 lags) spectrum (left) and lags vs frequency (right). Tiie tiiick solid lines represent the model 
REFLAGS for the best-fitting parameters (see Table 2, right column). 
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Figure 7. Obs. A (2004 spectrum (left) and lags vs frequency (right). The thick solid lines represent the model REFLAGS (including a 
Gaussian component in the case of the spectrum) for the best-fitting parameters (see Table 3, left column), in the frequency range 1-10 
Hz. 
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Figure 8. Obs. B' (2009 spectrum, 2004 lags) spectrum (left) and lags vs frequency (right). The thick solid lines represent the model 
REFLAGS for the best-fitting parameters (see Table 3, right column), in the frequency range 1-10 Hz. 



© 2012 RAS, MNRAS 000,[l}{TT] 



8 P. Cassatella et al. 



Parameter Obs. A + extra line Obs. B' 
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(9.11^=^5^) X 103 


Hout 1 Rout 


> 0.06 


< 0.05 


^src / Rg 


> 80 


> 80 


7 


> 1.1 


> 1.0 


i 


< 31.0 


< 30.2 


r 


1.459 ± 0.003 


1 4398+°-''°3 


log? 


< 1.34 


< 1.31 


Ec [keV] 


6.85 ±0.02 




<Tc [keV] 


0.25 ±0.03 




n==ff/(27r) 


0.81 


0.86 


(spectrum) 


1108 


1449 


(time-lags) 


74 


78 


xVdof 


1182 / 1402 


1527 / 1405 



Table 3. Spectral and lag fit parameters using the model REFLAGS 
plus an extra Gaussian (column one), and mixing the 2009 obser- 
vation spectrum with the 2004 observation lags in the range 1-10 
Hz. 




Energy [keV] 

Figure 10. Comparison between expected reprocessed emission 
using 2009 spectrum with 2004 lags (red) and 2009 lags (orange) 
and power-law fits to the data with photon indices —0.97 (same 
as low energy tail of reprocessed spectra, in blue) and 1.0 (as 
expected from synchrotron emission coming from a jet, in light 
blue) after removing extinction. The data points from UVOT are 
shown in green. 
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Figure 9. Power spectral densities for the 2004 (A, black) and 
2009 (B, red) observations, in the energy range 2.0 - 10.0 keV. 



that results (since the different line widths contributed by 
radii beyond a few thousand i?g cannot be resolved by the 
EPIC-pn detector). 

In a scenario where the lags are caused by reflection 
only at higher frequencies (e.g. with a propagation model 
explaining the lower- frequency lags) the lags would corre- 
spond to reflecting regions close to the central source of di- 
rect emission. If this is the case, excluding the lags below 
a certain frequency would lead to a lower and more plausi- 
ble i/out /flout- We therefore exclude lags below 1 Hz from 
the fit (Table 3) , since this frequency roughly corresponds 
to the threshold frequency below which the lags relative to 
softer energies are consistent with being due to propagation 
of disc fluctuations (Uttley et al. 2011). In addition, this 
range approximately corresponds to a range in frequency 
of the PSD (Fig. 9) where the PSD has a similar shape in 
both observations. Figs 7 and 8 show that when lags at fre- 
quencies < 1 Hz are excluded, the model cannot successfully 
reproduce either the lag shape or amplitude. 



4.1 Consistency of the reflection lags model with 
optical/UV data 

The parameters inferred for the best-fitting refiection model, 
which are quite extreme and still do not provide a good fit to 
the lags, can also be checked for the implied effect on optical 
and UV emission from GX 339-4. X-ray heating of the outer 
disc could in principle produce a large optical/UV flux if 
there is a large solid angle illuminated by the continuum, as 
is the case for the geometry inferred from our lag model fits. 
We can assume that each illuminated cell in the disc absorbs 
a fraction of the illuminating continuum equal to 1 — a{E), 
so that the incident luminosity absorbed by each cell can be 
calculated for the best-fitting given continuum shape and 
model ionisation parameter. If we make the simplifying as- 
sumption that the absorbed luminosity dominates over any 
intrinsic blackbody emission, we can equate the luminosity 
that is re-emitted by the cell to the absorbed luminosity and 
so determine the temperature of blackbody radiation emit- 
ted by each cell, and hence determine the total reprocessed 
contribution to the Spectral Energy Distribution (SED). 

To compare the predicted contribution to the opti- 
cal/UV SED from the geometry required by the lags model, 
we have extracted Swift/WOT (bands ubb, um2, uuu, 
uvv, uwl, uw2, uwh) spectra from a 1760 s observation of 
GX 339-4 made on 29 March 2009, two days after XMM- 
Newton observed the source. uvot2pha was used to extract 
spectra for source and background using a 6 arcsec radius, 
as well as extract response files. No additional aspect cor- 
rection was required. 

From the best-fitting spectral parameters found in fit B' 
and assuming a high-energy cut-off at 100 keV ( Motta et aQ 
2011), we derive an X-ray luminosity of 3.5 x 10'^' erg s~^ 
(assuming d — S kpc; Zdziarski et al. 2004). This value can 
be used to predict the reprocessed fiuxes that are consistent 
with the geometries inferred from the fits, accounting for 
interstellar extinction using the xspec model redden. 

Fig. 10 shows the expected reprocessed spectra for the 
two fits to the 2009 X-ray data (orange for B, red for B', 
black is the 2009 X-ray spectrum). The photon index F of 
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the power-law that characterises the reprocessed spectra at 
the energies covered by the UVOT data is -0.97. In these 
energy ranges, dust extinction needs to be taken into ac- 
count using the multiplicative model redden ( |Cardelli et al., 
1989[ >. E(B-V) = 0.933 is the value for the extinction cal- 
culated from IR dust maps along the line of sight towards 
our source^. By fitting a power-law with a photon index of 
-0.97 to the UVOT data, one finds the unabsorbed intrinsic 
power-law depicted in blue, which has a normalisation in the 
UVOT energy range that is several decades larger than that 
expected from reprocessing, and requires E(B-V) = 1.587 
for x^/dof=0.8/4. Therefore the model severely underpre- 
dicts the observed flux. A more likely explanation for the 
optical/UV emission is flat-spectrum synchrotron emission 
from a compact jet ||Maitra et al.ll2009|, or magnetised hot 
accretion flow (Vclcdina et al. 2011). Assuming a power-law 
photon index of 1.0 (spectral energy index of 0), one also 
obtains a good fit (x^/dof=0.64/4) and a lower extinction 
than in the previous case, E(B-V) = 1.089 (light blue line). 
This is consistent with the results found by [Maitra et al.| 
(2009) for the same object, fitting broadband data using 
only synchrotron and inverse-Compton models. Therefore, 
the UVOT data cannot be explained solely by reprocessing 
in the flared disc envisaged by the lags model and is more 
likely to be produced by a synchrotron process. However, a 
small contribution from reprocessing cannot be ruled out. 



5 DISCUSSION 

5.1 Physical implications 

The analysis shown in Section 4 shows that extreme disc 
geometries arc preferred for the model to reproduce the ob- 
served lags as closely as possible within the constraints given 
by the overall spectral shape, including the strength of the 
iron feature at 6.4 keV. There are several effects of the re- 
processing geometry on the spectral and lag data that need 
to be highlighted to understand why the spectrum can be 
fitted well, whereas the lags cannot. 

Firstly, while the amount of flux in the iron line is pro- 
portional to the solid angle subtended by the disc, the abil- 
ity to determine how much of it is produced in the outer 
radii of the disc (where Doppler effects are weak) is limited 
by the resolution of the XMM-Newton EPIC-pn detector. 
Therefore, the description of the geometry that could be in- 
ferred by the spectral modelling alone is degenerate, since 
line emission from the largest radii (e.g. ~ 10^ -Rg) cannot 
be resolved from that at more modest (but still large) radii 
(e.g. ~ 10'' Rg). The result of this effect is that the spectral 
fits are not sensitive to variations in the outer radius of the 
reflector at large radii. 

On the other hand, the lags are very sensitive to the 
geometry at large radii. Firstly, the lags at low frequencies 
increase with both the solid angle and light-travel time to 
the reflector at large radii. A larger outer radius and more 
flared disc therefore corresponds to larger lags. However, the 
size-scale of the largest radius also corresponds to a charac- 
teristic low-frequency flattening in the lag versus frequency 
dependence. This is because the frequency-dependent drop 
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in lags seen at higher frequencies is caused by smearing of the 
reflection variability on time-scales shorter than the light- 
travel size-scale of the reflector. The reflection variability 
amplitude is not smeared out for variability time-scales sig- 
nificantly longer than the light-travel time to the largest 
disc radii, and the lag at low frequencies quickly approaches 
the average light-travel delay from the reflector (diluted by 
the direct continuum emission which has zero intrinsic lag). 
The frequency of this characteristic flattening in the lag- 
frequency dependence is therefore a sensitive indicator of the 
size-scale of the reflector. However, it is not possible to rec- 
oncile the position of this flattening at frequencies ~ 0.2 Hz 
with the large amplitude of the lags at low frequencies, which 
imply an even larger reflector subtending an even greater 
solid angle to the continuum at large radii. 

The result of those effects is that despite the already 
extreme inferred geometries the model is clearly underprc- 
dicting the lags. At this point, the maximum value of the 
lags is now constrained by the spectral modelling, which 
cannot place a tight constraint on the disc outer radius but 
does limit the solid angle of the reflector. It is instructive 
to consider the effects on the predicted spectral shape when 
the model is fitted to the 2004 lags alone (fixing F = 1.5) 
and the same parameters are used to estimate the resulting 
spectrum. This yields x'^/<iof=66. 4/16=4. 15 for the fit to 
the lags and yields an apparent solid angle fl'^'^ /2tt — 1.51 
subtended by the disc. The resulting spectral shape is com- 
pared to the data in Fig. [ill Therefore, even fitting the lag 
data alone with the model cannot produce a good fit, and 
the fit that is obtained shows that much larger solid angles 
of large-scale refiection are required than are permitted by 
the spectrum. 

The lags cannot be explained solely by reflection, it is 
therefore necessary to invoke an additional mechanism to 
explain them. This result is perhaps not surprising, since we 
have previously found evidence that in hard state BHXRBs, 
fluctuations in the accretion disc blackbody emission are cor- 
related with and precede the variations in power-law emis- 
sion (Uttley et al. 2011). Although the disc variations seem 
to drive the power-law variations, this does not in itself ex- 
plain the lags within the power-law band, which we consider 
in this paper, since the disc emission only extends up to 
~ 2 keV. However, as we noted in 'Uttley et al. (2011), at 
frequencies < 1 Hz, the lags of the power-law emission rela- 
tive to the disc-dominated 0.5-0.9 keV band show a similar 
frequency-dependence to the lags seen within the power-law 
band (i.e. r oc i/^"'^). This strongly suggests that the lags 
intrinsic to the power-law are somehow connected to the 
mechanism which causes the power-law to lag the disc, most 
probably due to the propagation of accretion fluctuations 
through the disc before reaching the power-law emitting hot 
flow. 

One possibility is that the disc is sandwiched by the 
hot-flow/corona which produces power-law emission which 
becomes harder towards smaller radii, leading to hard lags 
Eis fluctuations propagate inwards (e.g. iKotov et al.||2001| 
Arcvalo & Uttley 2006). This model can explain the hard 
lags in terms of propagation times in the flow, which are 
much larger than light-crossing times and so can produce 
relatively large lags which the refiection model struggles to 
produce without leading to solid angles of large-scale refiec- 
tion. Reflection may also contribute to the lags at some level, 
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but is not the dominant mechanism, at least at frequencies 
< 1 Hz. 

A more detailed analysis of the contribution of reflec- 
tion to the observed lags could be performed using datasets 
with higher signal-to-noise, by e.g. searching for reflection 
signatures around the iron line investigated by Kotov et al. 
(2001). Lags vs. energy spectra for GX 339-4 are shown in 
[Uttley et aL\ |2011| ) for the 2004 XMM-Newton observation, 
and demonstrate that the current quality of the data is not 
sufficient to detect these features. 

It is also possible that multiple distinct components 
contribute to the lag, e.g. associated with the different 
Lorentzian features that contribute to the PSD. This pos- 
sibility could explain the apparent 'stepping' of the lag vs. 
frequency that appears to be linked to the frequencies where 
the dominant contribution to the PSD changes from one 
Lorentzian component to another (Nowak 2000 ). 
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Figure 11. Spectral model (reflags) compared to 2004 spec- 
trum using the parameters that are derived from fitting the 2004 
lags alone. The lower panel shows the ratio data/model. 



5.2 Wider implications of combined 
spectral-timing models 

In this work we have considered (and ruled out) a relatively 
simple model for the lags in terms of the light-travel times 
from a large scale reflector. However, it is important to stress 
the generalisability of our approach to other models. In par- 
ticular, we have shown how it is possible to combine timing 
and spectral information to fit models for the geometry and 
spatial scale of the emitting regions of compact objects. Pre- 
vious approaches to use the information from time-lags to fit 
models for the emitting region have focussed on fitting the 
lag data (e.g. Kotov et al. 2001; Poutanen 2002). However, 
since these models also make predictions for the spectral be- 
haviour, it is possible to achieve stronger constraints on the 
models by fitting the lags together with the time-averaged 
energy spectrum, as done here, or with spectral-variability 
products such as the frequency-resolved rms and covariance 
spectra ( Revnivts ev et al.ylQQQi [Wilkinson fc Uttleyn2009( 
jUttley et al. 2011). 

In order to use these techniques more generally, one 
needs to calculate the energy-dependent response function 
for the emitting region, i.e. determine the emission as a func- 
tion of time delay and energy. This approach can be used to 
test reverberation models for the small soft lags seen at high 
frequencies in AGN (Fabian et al. 2009, Zoghbi et al. 2010; 
de Marco et al. 2011; Emmanoulopoulos et al. 2011) and 
BHXRBs (Uttley et al. 2011), which offers the potential to 
map the emitting region on scales within a few gravitational 
radii of the black hole. Other models can also be considered, 
e.g. to test the propagation models for the low-frequency lags 
with the time-delay expressed in terms of propagation time 
through the accretion fiow. Future, large-area X-ray detec- 
tors with high time and energy-resolution, such as the pro- 
posed ATHENA and LOFT missions, will allow much more 
precise measurements of the lags in combination with good 
spectral measurements, so that fitting of combined models 
for spectral and timing data could become a default ap- 
proach for studying the innermost regions of compact ob- 
jects. Future research in this direction is strongly encour- 
aged. 
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